Theory of pairing symmetry inside the Abrikosov vortex core 
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We show that the Cooper pair wave function at the center of an Abrikosov vortex with vorticity 
m has different parity with respect to frequency from that in the bulk if m is an odd number and has 
the same parity if m is an even number. As a result, in a conventional vortex with m = 1, the local 
density of states at the Fermi energy has a maximum (minimum) at the center of the vortex core 
in even(odd)-frequency superconductor. We propose a scanning tunneling microscope experiment 
using a superconducting tip to explore odd-frequency superconductivity. 



The study of the mixed state in type-II superconduc- 
tors has a long history and revealed a variety of physi- 
cal phenomena. [l[ In the clean limit, low-energy bound 
states (the Andreev bound states) are generated in the 
vortex core due to the spatial structure of the supercon- 
ducting pair potential. [1, Q One of the manifestations 
of the bound states is the enhancement of zero-energy 
quasiparticle density of states (DOS) locally in the core, 
observable as a zero-bias conductance peak by scanning 
tunneling microscope (STM).[1,I3 However, despite ex- 
tensive studies of the vortex core, the issue of pairing 
symmetry in the core remains unexplored. 

Generally, superconducting pairing is classified into 
even-frequency or odd-frequency state according to a 
symmetry with respect to time. Due to the Fermi 
statistics, even-frequency superconductors belong to the 
symmetry class of spin-singlet even-parity (ESE) or 
spin-triplet odd-parity (ETO) pairing state, while odd- 
frequency superconductors belong to the spin-singlet 
odd-parity (OSO) or spin-triplet even-parity (OTE) pair- 
ing state. Although vortex core state of even-frequency 
superconductors has been well studied, that in odd- 
frequency superconductors has not been clarified yet. 

The possibility of the odd-frequency pairing state in 
various kinds of uniform systems was discussed in Refs. 
|3, albeit its realization in bulk materials is still 
controversial. On the other hand, the realization of 
the odd-frequency pairing state in inhomogeneous even- 
frequency superconducting systems has recently been 
proposed. It is established that odd-frequency pairing 
is induced due to symmetry breaking in such systems. 
In ferromagnet/superconductor junctions, odd-frequency 
pairing emerges due to the broken symmetry in a spin 
space. [7| It was recently realized that proximity- induced 
odd-frequency pairing may be generated near normal 
metal/superconductor interfaces due to the breakdown of 
translational symmetry 0] or in a diffusive normal metal 
attached to a spin-triplet superconductor. [9( 

Since an Abrikosov vortex breaks translational sym- 
metry in a superconductor, one may expect the emer- 
gence of an odd-frequency pairing state around the vor- 
tex core even in a conventional spin-singlet s-wave su- 
perconductor. On the other hand, one may imagine that 



even-frequency pairing state may be generated around 
the vortex core in odd-frequency superconductor. 

In this Letter, based on the quasiclassical theory of su- 
perconductivity, we develop a general theory of pairing 
symmetry in an Abrikosov vortex core in clean supercon- 
ductors, including odd-frequency superconductivity. We 
show that for a vortex with vorticity m in superconduc- 
tor, pairing function of the Cooper pair at the vortex 
center has the opposite (same) symmetry with respect to 
frequency to (as) that of the bulk if m is an odd (even) 
integer. For a conventional vortex with m = 1, we show 
that zero energy local DOS is enhanced (suppressed) at 
the center of the vortex core in even (odd)-frequency su- 
perconductor. We further reveal that OSO p-wave pair- 
ing is generated at the center of the core of ESE s-wave 
superconductor. On the other hand, in OSO p-wave su- 
perconductor, ESE s-wave pairing state emerges at the 
center of the vortex core. Based on these results, we 
propose an experimental setup to explore odd-frequency 
superconductivity by probing a local Josephson coupling 
by STM with superconducting tip. 

The electronic structure of the vortex core in a single 
Abrikosov vortex in a clean superconductor is described 
by the quasiclassical Eilenberger equations 13, U | based 
on the Riccati parametrization of the quasiclassical prop- 
agator Along a trajectory r(x') = r + x' vp with 
unit vector vp parallel to vp, the Eilenberger equations 
are generally represented in 4x4 matrix form.[13j| For a 
singlet (triplet) superconductor with A = Aa y (a x ) [a x 
and <jy are Pauli's matrices in spin space) [lij . these equa- 
tions are reduced to the set of two decoupled differential 
equations of the Riccati type for the functions a{x') and 
b(x'), 

hv F d x ,a{x') + [2e n + A^a{x')] a(x') - A = 0, 
hv F d x ,b{x') - [2e n + Ab(x')} b(x') + A f = (1) 

where ie n are the Matsubara frequencies, and A T = 
(— )A* for even (odd)-frequency superconductor. For 
simple case of a cylindrical Fermi surface, the Fermi ve- 
locity can be written as vp = vp(ei cosO + e 2 sin.0). 
We choose the following form of the pair potential: 



A(r, 6, E) = A o x(0, E)F{r) exp(im<p) 



(2) 
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with r — x 2 + y 2 and exp(itp) = (x + iy)/ ' \J x 2 + y 2 . 
Here, F(r) denotes the spatial profile of the gap, to is the 
vorticity, and E) is the symmetry function. Also, we 
introduce the coherence length £ = Tlvf/Aq, the center 
of a vortex is situated at x — y — 0, and exp(irrup) is the 
phase factor which originates from the vortex. 

We obtain pairing function of the Cooper pair (anoma- 
lous Green's function) / as / = — 2a/ (1 + ab). For the 
calculation of the local DOS normalized by its value in 
the normal state, the quasiclassical propagator has to be 
integrated over the angle 9 which defines the direction of 
the Fermi velocity. The normalized local DOS in terms 
of functions a and b is given by 



TABLE I: Pairing symmetry in the vortex state. 



N(r ,E) 





'1 - ab' 




1 + ab 


/O 27T 



(3) 
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where E denotes the quasiparticle energy with respect 
to the Fermi level and <5 is an effective scattering param- 
eter that corresponds to an inverse mean free path. In 
numerical calculations throughout this paper, we will fix 
this value as 5 = 0.1A . 

First, we discuss the general property of the symme- 
try at the vortex center. Vorticity and the symmetry 
of superconductor with respect to frequency crucially af- 
fect the symmetry of the Cooper pair at the core center. 
Consider a trajectory passing through the center of the 
vortex. By setting x' = at the vortex center, we get 
b{x',t n ) = —l/a(—x',—e n ) from the Eilenberger equa- 
tions for even-frequency superconductor with odd integer 
to or odd-frequency superconductor with even integer to. 
Similarly, we obtain b{x',t n ) — l/a(—x',—e n ) for odd- 
frequency superconductor with odd integer m or even- 
frequency superconductor with even integer to. Thus, 
at the vortex center x 1 = 0, we get f(e n ) — — /(— e„) 
in the former case, while /(e„) = /(— e n ) in the latter. 
Note that spin is conserved in the vortex state consid- 
ered. Therefore, quite generally, for an odd integer to the 
induced pairing at the vortex center has different sym- 
metry with respect to frequency from that in the bulk 
superconductor. On the other hand, for an even integer 
m, the induced pairing at the vortex center has the same 
symmetry as that of the bulk. We summarize pairing 
symmetry at the vortex center in Table I. 

For conventional s-wave case, there have been several 
studies of multi- vortex state with to > 1 , It was 
shown that zero energy peak in the DOS only appears for 
odd number to at the vortex center [l6| • This statement 
is consistent with our result for the conventional s-wave 
case of x(6*, E) = 1 because odd-frequency pairing state is 
generated only for odd integer to. The relation between 
zero energy peak in DOS and odd-frequency pairing state 
will be discussed later. 





bulk state 


vorticity 
m 


symmetry at 
the center 


(1) 


ESE 


odd 


OSO 


(2) 


ESE 


even 


ESE 


(3) 


ETO 


odd 


OTE 


(4) 


ETO 


even 


ETO 


(5) 


OSO 


odd 


ESE 


(6) 


OSO 


even 


OSO 


(7) 


OTE 


odd 


ETO 


(8) 


OTE 


even 


OTE 



0) 
T3 

Q. 

E 

CO 
s_ 

CD 
D_ 




0.4 0.6 0. 



FIG. 1: (color online) Results for ESE s-wave superconduc- 
tor, (a) Normalized local DOS around the vortex at E — 0. 
The center of the vortex is situated at x = y = 0. Spatial 
dependencies of (b) ESE and (c) OSO components at E = 0. 
Only OSO components, Kef px and Im/ ra , can survive inside 
the core (near x — 0) . 



In general, the most realizable vorticity is m = 1. 
Thus, in the following, we will study in detail two typ- 
ical cases at to = 1 with ESE s-wave and OSO p x - 
wave superconductors where we choose x(#, E) = 1 and 
X{6, E) = (C cos 9E/A )/[l + {E/A ) 2 } with C = 0.8,0 
respectively. Also, spatial dependence of the gap is cho- 
sen as F(r) = tanh(r/£). 
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FIG. 2: (color online) Results for OSO p^-wave superconduc- 
tor, (a) Normalized local DOS around the vortex at E — 0. 
The center of the vortex is situated at x — y = 0. Spatial 
dependencies of (b) OSO and (c) ESE components at E — 0. 
Only ESE components can exist inside the core (near x = 0). 



Due to the broken translational symmetry of the sys- 
tem, various pairing states are expected to emerge around 
the vortex. In order to study possible pairing states, we 
decompose anomalous Green's function / into various an- 
gular momentum components as follows 



/ - fa + fpx cos 9 + f py sin 9 + f dx 2_ y 2 cos 29 
+fdx V sin 29 + ffi cos 39 + // 2 sin 39 + 



(4) 



Note that all the above pairing components are spin sin- 
glet. 

Figure [1] shows the results for ESE s-wave supercon- 
ductor. The local DOS around the vortex at E = 
is shown in Fig. [T] (a). As is well known, zero energy 
peak appears in the core Q. The spatial dependencies 
of decomposed anomalous Green's function / at E = 
are shown in Figs. [2(b) and (c). Interestingly, only OSO 
pairing components, Kef px and lmf py , survive at the cen- 
ter of the core. With the increase of the distance from 



the core center, the magnitudes of / decrease rapidly, 
except for s-wave one. Note that other angular momen- 
tum components not shown in this figure are negligibly 
small. We also find that the following representation of 
anomalous Green's function / holds at the center of the 
core: 

/ = (Rcf px + ibnfpx) (cos 9 + i sin 9) = f px exp(£0). (5) 

Thus, we see that anomalous Green's function at the core 
center has chiral p-wave symmetry. 

The enhancement of the local DOS in the presence 
of odd-frequency pairing can be understood, irrespective 
of the detailed shape of A by using the normalization 
condition for the quasiclassical Green's functions, g 2 + 
f f = 1. Indeed, since for odd- frequency pairing state, 
the anomalous Green's function / = — 26/(1 + ab) at E = 
is given by J (9) = -f*{9) (see Ref.[I3) and local DOS 
is given by N(E) = —Keg, one can show that generally 
N(E = 0) > 1 since g 2 = 1 + |/| 2 > 1. This means that 
the emergence of the odd-frequency pairing is a physical 
reason of zero energy peak of the local DOS inside the 
core. The manifestation of the odd-frequency chiral p- 
wave pairing state at the center of the vortex core is also 
consistent with the experimental fact that the observed 
zero-bias conductance peak by STM at a vortex center 
is very sensitive to disorder [171 ]. since p-wave pairing is 
fragile to impurity scattering. 

Figure [2] depicts the results for OSO p^-wave super- 
conductor. The local DOS around the vortex at E = is 
shown in Fig. [5] (a) . In dramatic contrast to the result for 
the ESE s-wave, zero energy DOS is suppressed at the 
core. The spatial dependencies of decomposed anoma- 
lous Green's function at E = are shown in Figs. [2] (b) 
and (c). As is seen, only ESE pairing components ex- 
ist at the center of the core. For even-frequency pairing 
state, f(9) — f*{9) is satisfied at E = and hence we 
get N(E = 0) < 1 which is consistent with Fig. [2] (a). 
By comparing Figs.l and 2, it is clear that zero energy 
local DOS N(0) has a maximum at the center of the vor- 
tex core in even-frequency superconductor, while it has a 
minimum at the core center in odd-frequency supercon- 
ductor. This difference can be detected by STM. 

As regards the candidate for the odd-frequency super- 
conductor, CeCu2Si2 and CeRhlns are possible materials 
0, Eij . In these systems, OSO state with p-wave sym- 
metry is considered to be promising Q. In the light of 
the present theory, ESE s-wave pairing is expected to 
appear inside the vortex core. Based on this idea, we 
propose an experimental setup to verify the existence of 
odd-frequency pairing in bulk materials by using super- 
conducting STM where we use conventional s-wave su- 
perconductor as a STM tippjl as shown in Fig. 3. Local 
Josephson current measured in STM experiment with su- 
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FIG. 3: (color online) Suggested experimental setup to probe 
the local Josephson coupling between superconducting tip and 
superconductor with Abrikosov vortex core. 

perconducting tip is given by (20j 

eIR = ttTJ2 Im (/*(#, e n )f s (9, £„))• (6) 

Here, R is the junction resistance, T is temperature, and 
fs is anomalous Green's function in the STM tip. It fol- 
lows from this expression that finite Josephson current is 
allowed only when superconducting STM has the same 
symmetry as that in the vortex state. Therefore, finite 
Josephson current is allowed at the vortex core only when 
the pairing symmetry at the core is ESE s-wave. First, we 
consider bulk superconductor which has an ESE s-wave 
symmetry with m=l. As seen from Fig. [TJ only OSO 
state is generated at the core center. Then, the local 
Josephson coupling is absent at the core while it exists 
in the bulk. On the other hand, in OSO p-wave supercon- 
ductor with m = 1 , ESE s-wave state is generated inside 
the core (see Fig. In this case, the local Josephson 
coupling only exists at the core. Such a dramatic feature 
is useful to explore and identify odd-frequency pairing in 
bulk superconductor. We summarize the possible local 
Josephson couplings with ESE s-wave STM tip in Table 
II. 



TABLE II: Local Josephson coupling with ESE s-wave STM 
tip. 



Bulk state 


Local Josephson coupling 


core 


bulk 


ESE 


No 


Yes 


OSO 


Yes 


No 



In summary, we have developed a general theory of 
pairing symmetry inside the Abrikosov vortex core in su- 
perconductors, including odd- frequency superconductiv- 
ity. We have found that for a vortex with vorticity m 
in a superconductor, anomalous Green's function at the 
vortex center has the opposite (same) symmetry with re- 
spect to frequency to (as) that of the bulk if m is an odd 



(even) integer. We have also shown that zero energy local 
density of states is enhanced (suppressed) at the center of 
the vortex core for even(odd)-frequency superconductor. 
Based on the obtained results, we proposed a scanning 
tunneling microscope experiment using a superconduct- 
ing tip to detect local Josephson coupling in order to 
explore and identify odd-frequency superconductor. 
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